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ABSTRACT 

I present a XMM-Newton observation of the accretion driven millisecond X-ray pulsar SAX 
J1808. 4-3658 during its 2000 outburst. The source was conclusively detected, albeit at a level of 
only ^ 2 X lO'^'^ erg s~^. The source spectrum could be fitted with a power-law model (with a 
photon index of ^^2.2), a neutron star atmosphere model (with a temperature of ~0.2 keV), or 
with a combination of a thermal (either a black-body or an atmosphere model) and a power-law 
component. During a XMM-Newton observation taken approximately one year later, the source 
was in quiescence and its luminosity was a factor of ^4 lower. It is possible that the source 
spectrum during the 2000 outburst was softer than its quiescent 2001 spectrum, however, the 
statistics of the data do not allow to make a firm conclusion. The results obtained are discussed 
in the context of the 2000 outburst of SAX J1808. 4-3658 and the quiescent properties of the 
source. 

Subject headings: accretion, accretion disks — stars: individual (SAX J1808. 4-3658) — stars: neutron 
— X-rays: stars 



1. Introduction 

The X-ray transient SAX J1808.4-3658 was dis- 
covered in September 1996 when it exhibited a 
weak outburst lasting only a few weeks (In 't Zand 
et al. 1998, 2001). In April 1998 the source was 
found to be in outburst again (Marshall 1998) and 
it was discovered that the source exhibits coherent 
millisecond X-ray oscillations with a frequency of 
approximately 401 Hz (Wijnands & van der Klis 
1998). In early 2000, the source exhibited a third 
outburst during which it showed erratic luminos- 
ity behavior with luminosity swings of three or- 
ders of magnitude within a few days (Wijnands et 
al. 2001, 2002). This erratic behavior lasted for 
several months before the source returned to qui- 
escence. Very recently, in October 2002, a fourth 
outburst of the source was detected (Markwardt, 
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Miller, Wijnands 2002) during which its peak lu- 
minosity was very similar to that observed during 
the 1996 and 1998 outbursts. 

In quiescence, SAX J1808. 4-3658 has been ob- 
served on several occasions with the BeppoSAX 
and ASCA satellites (Stella et al. 2000; Dotani, 
Asai, & Wijnands 2000; Wijnands et al. 2002). 
The source was very dim in quiescence, with a lu- 
minosity close to or lower than 10"^^ erg s^^. Due 
to the limited angular resolution of BeppoSAX, 
doubts were raised as to whether the source de- 
tected by this satellite was truly SAX J1808.4- 
3658 or an unrelated field source (Wijnands et al. 
2002). Campana et al. (2002) reported on a qui- 
escent observation of the source performed with 
XMM-Newton which resolved this issue. They de- 
tected the source at a luminosity of 5 x 10^^ erg 
s-i and found that the field around SAX J1808.4- 
3658 is rather crowded with weak sources. Two 
such sources are relatively close to SAX J1808.4- 
3658 and might have conceivably caused a system- 
atic positional offset during the BeppoSAX obser- 
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vations of SAX J1808.4 3658. 

After it was found that the source had become 
active again in January 2000 (Van der Klis et 
al. 2000), a Director's Discretionary Time (DDT) 
observation request was submitted to the direc- 
tor of XMM-Newton to study the outburst X-ray 
spectrum of SAX J1808.4 3658. This request was 
granted and on March 6, 2000, XMM-Newton ob- 
served the source. Due to the extreme variabil- 
ity of the source during its 2000 outburst (Wij- 
nands et al. 2001), the XMM-Newton observation 
was performed during times when the source had 
very low luminosities (see § 2). Because of this 
the DDT observation was considered to be in con- 
flict with an approved Cycle 1 observation on this 
source with the purpose of studying the quiescent 
X-ray properties of SAX J1808.4-3658 (see Cam- 
pana et al. 2002, for the results of that observa- 
tion). Consequently, the DDT observation was not 
made public until recently. In this paper I discuss 
this DDT observation during which the source was 
conclusively detected at a luminosity of ^ 2 x 10^^ 
erg s~^. 

2. Observation and analysis 

Because a bright X-ray source was expected 

during the 2000 outburst, the XMM-Newton/F.PlC- 
pn camera was used in the timing mode to elimi- 
nate pile-up problems and to study the pulsations. 
However, due to the faintncss of the source and 
the high background when using this mode, the 
resulting data are not optimal for studying the 
source spectrum at the observed low luminosities. 
Therefore, only the data obtained with the two 
metal oxide semiconductor (MOS) cameras will 
be discussed. The MOSl camera was used in the 
Full Frame mode, but the M0S2 was in the Small 
Window mode (again to limit the anticipated pile- 
up). For both cameras the medium filter was used. 
The source was not detected in the RGS instru- 
ments so that data will not be discussed further. 
The log of the observation and the data used in 
this paper are listed in Table 1. 

The data were analysed using the Standard 
Analysis System (SAS), version 5.3. The ob- 
servation was split into two observation identi- 
flers (IDs; see Tab. 1) and it was found that the 
MOSl science data were all present in the first 
ID (0119940201), but that the housekeeping data 
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Fig. 1. — XMM-Newton/MOSl images of the re- 
gion near SAX J1808.4-3658, which is the bright- 
est source in the images. North is up and East 
is to the left. The BeppoSAX, ASCA/GIS, and 
^5'Cj4/SIS Ifj error circles of the sources detected 
by those instruments are shown (Wijnands et al. 
2002). The energy band corresponding to each 
image is shown on the bottom left. 



were split over the two IDs (the M0S2 science data 
anrf housekeeping data were both split, in a consis- 
tent manner, over the two IDs). This discrepancy 
caused the standard processing of the data to cre- 
ate only about half the available MOSl science 
data. To obtain the full amount of MOSl science 
data, the SAS task emchain was executed but after 
the housekeeping files for the MOSl were merged^ 
and using the latest calibration files. Those newly 
created event list files (also for the M0S2) were 
used in the subsequent analysis. At the end of the 
observation, a strong background fiare occurred. 
The data in which the total count rate (all CCDs 
and no extra filtering) exceeded 5 and 4 counts per 
second (using 100 second bins) for the MOSl and 

^ Using the threads listed at 

http://wave.xray.mpe.mpg.de/xmm/cookbook/EPIC_PN/merge_odf.html 
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M0S2 data, respectively, were excluded from the 
analysis. The resulting total live time of the cen- 
tral CCD (the CCD on which SAX J1808.4-3658 
is located) is listed in Table 1. 

The MOSl image'' in the 0.5-10 keV energy 
band is shown in Figure 1 (top) where we clearly 
see that SAX J1808.4-3658 was detected. This 
source is the brightest source on the central CCD, 
strongly indicating that during the BeppoSAX 
2000 outburst observations most of the detected 
flux originated from SAX J1808. 4-365 and not 
from an unrelated field source as was suggested 
by the BeppoSAX data (Wijnands et al. 2002). 
Campana et al. (2002) reached similar conclusions 
using a quiescent observation of the source and 
they suggested that the systematic offset of the 
BeppoSAX position with regards to that of SAX 
J1808. 4-3658 might have been caused by the pres- 
ence of two faint sources close to SAX J1808.4- 
3658 (Fig. 1, top panel; see also Fig. 1 in Campana 
et al. 2002). Moreover, due to the source faint- 
ness during the BeppoSAX observations, only the 
Medium Energy Concentrator Spectrometer pro- 
duced useful data (Wijnands et al. 2002) and this 
instrument was only sensitive in the energy range 
1.3-10 keV. Therefore, images of the data in the 
0.5-1.3 keV and the 1.3-10 keV energy bands were 
made (Figs. 1, middle and bottom panel) and SAX 
J1808. 4-3658 is detected in both energy ranges. 
In contrast, the two extra sources are only de- 
tected in the 1.3-10 keV band, demonstrating that 
the fractional flux contribution of the two other 
sources to the combined flux increases with pho- 
ton energy. These differences in source spectra add 
further evidence to the suggestion that those two 
extra sources might have caused the systematic 
offset in the BeppoSAX observations. It should be 
noted that the BeppoSAX fluxes quoted by Wij- 
nands et al. (2002) and Stella et al. (2000) for SAX 
J1808.4 3658 arc likely close to its true flux. A 
small contamination of the fluxes by the two close- 
by sources is likely, but those sources are consid- 
erably loss luminous than SAX .11808.4-3658 and 
therefore the contamination should be small. 

The X-ray source spectra were extracted using 
a circle with a radius of 20" on the position of 

*The use of the Small Window mode for the MOS2 camera 
limited the size of the image around SAX J 1808.4-3658 and 
therefore this image is not displayed. 
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Fig. 2.— XMM-Newton/MOSl (circles) and 
M0S2 (triangles) spectra of SAX J1808.4-3658 
during the 2000 outburst. The solid lines repre- 
sent the best power-law fit to the data. 



SAX J1808.4-3658. For the MOSl camera, the 

background spectrum was extracted using a cir- 
cle with a radius of 200" on the same position, 
but excluding the detected point sources in this 
region. Because of the use of the Small Window 
mode, only a limited field could be use to extract 
the background spectrum for the M0S2 camera 
(i.e., an annulus was used on the source position 
with an inner radius of 30" and an outer radius of 
50")^. The RMF and ARF files were created with 
the SAS tools rmfgen and arfgen. The spectra ob- 
tained were then grouped using the FTOOL grp- 
pha into bins with a minimum of 20 counts per bin 
to validate the use of the statistics. The MOSl 
and M0S2 spectra are shown in Figure 2 and 
were fitted (using XSPEC version 11.1.0; Arnaud 
1996) simultaneously using the same model (sec 
Tab. 2 for the fit parameters). The column density 
A^H was allowed to float and the value obtained 
was always consistent with the value (1.22 x lO^"' 
cm~^) inferred from the Ay measured by Wang et 

^The accuracy of the background subtraction was checked 
by using background regions located on the other, fully 
read-out, CCDs. The fluxes were always within 6% of each 
other regardless of the background used, and the spectral 
parameters within 2%. To avoid systematic uncertainties 
due to the different responses and different offsets of the 
individual CCDs, the background obtained from the central 
CCD was preferred. 
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al. (2001) (and using the relation between TVh and 
Ay from Predehl & Schmitt 1995). The spectrum 
could be fitted with a power-law model with in- 
dex of 2. 2 ±0.3 and a 0.5 10 keV luminosity of 
1.7 X 10^^ erg s~^ (assuming a distance of 2.5 
kpc; In 't Zand et al. 2001). A black body model 
could not fit the data accurately (x^ /degrees of 
freedom = 88.7/25). A neutron star atmosphere 
model (that of Zavlin, Pavlov, & Shibanov 1996) 
could fit the data but with a relatively high tem- 
perature kT of ^^0.2 keV (for an observer at infin- 
ity) and a neutron star radius of 4.5±0.1 km (as 
measured on the surface and using a neutron star 
mass of 1.4 Mq). 

The X-ray spectra of other quiescent neutron 
star transients can sometimes be described by 
a two component model (a soft thermal compo- 
nent below 1 keV and a power-law-like component 
above a few keV; e.g., Asai et al. 1996, 1998) but 
with the soft component dominating the spectrum 
(although the power-law component can occasion- 
ally contribute nearly half the 0.5-10 keV flux of 
the source). Although not required by the data, 
the spectra were also fitted using such a two com- 
ponent model, with either a black body or an at- 
mosphere model for the soft component. With 
these models, the spectra could be accurately fit- 
ted (see Tab. 2), although the atmosphere plus 
power law combination was very unstable and the 
errors on the parameters are therefore relatively 
large. The temperatures obtained for the soft com- 
ponent were similar to those obtained for other 
quiescent neutron star transients. The flux con- 
tribution of the soft component to the 0.5-10 keV 
flux was only ~25% of the total flux in contrast 
to the other systems in which the soft component 
dominates. 

3. Discussion 

This paper discusses the XMM-Newton obser- 
vation of the accretion-driven millisecond X-ray 
pulsar SAX J1808. 4-3658 performed during its 
2000 outburst. Similar to the BeppoSAX obser- 
vations performed around the same time (one of 
those observations was on the same day as the 
XMM-Newton one), the XMM-Newton observa- 
tion revealed only a weak X-ray source with a 0.5- 
10 keV Imninosity of ^ 2 x lO"*^ erg at the posi- 
tion of SAX J1808. 4-3658. Despite its weakness. 



the source was the brightest one on the central 
CCD proving that the BeppoSAX source detected 
by Stella et al. (2000) and Wijnands et al. (2002) 
is indeed SAX .J1808.4 3658 and not an unrelated 
field source. The systematic off-set between the 
measured and the true position of SAX J1808.4- 
3658 in those BeppoSAX observations can likely 
be explained by two nearby field sources. These 
sources might have also contaminated the flux as- 
signed to SAX J1808. 4-3658 (see also Campana et 
al. 2002). 

Campana et al. (2002) reported on a quies- 
cent observation of SAX J1808. 4-3658 performed 
with XMM-Newton about a year after the end of 
the 2000 outburst (the spectral fit parameters re- 
ported by Campana et al. (2002) are also listed 
in Tab. 2 for comparison). They found that the 
source had a luminosity of only ~ 5 x 10^^ erg 
s"-"-, about a factor of 4 lower than what was 
measured during the 2000 outburst observation. 
This demonstrates that at very low luminosities, 
SAX J1808.4 3658 can exhibit variability and in- 
dicates that the source probably displays variabil- 
ity in its quiescent state. However, caution is ad- 
vised because the 2000 outburst observation was 
performed when the source exhibited violent be- 
havior (Wijnands et al. 2001). Therefore, it is 
still possible that in "true" quiescence the source 
will be observed consistently at the low level re- 
ported by Campana et al. (2002). The spectrum 
of SAX J1808.4-3658 during the 2000 outburst ap- 
pears softer than its quiescent spectrum: when 
the spectrum is fitted with a power-law model, 
the photon index is 2.2±0.3 during the 2000 out- 
burst vs. 1.5lo!3 in quiescence. The fact that 
the 2000 outburst spectrum can be accurately fit 
with a neutron star atmosphere model or with 
a combination of a thermal plus power-law com- 
ponent (in contrast with the quiescent spectrum; 
Campana et al. 2002), also suggests a difference 
in the source spectra between the two epochs. It 
is interesting to note that in the two-component 
model, the photon index during the 2000 outburst 
observation is very similar to that measured dur- 
ing the quiescent observation, suggesting that the 
shape of the power-law component might not have 
changed considerably between the different epochs 
(although the flux of this component was still a 
factor of ^3 higher during the 2000 outburst than 
in quiescence) . Due to the limited statistics of the 
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data, the exact degree of spectral variability can- 
not be determined nor can the exact cause of these 
possible variations be established. 

During the 2000 outburst, the source fluctuated 
in luminosity by over 3 orders of magnitude on 
timescales of days (Wijnands et al. 2001). How- 
ever, the BeppoSAX observations reported by Wi- 
jnands et al. (2002) provided only rough estimates 
of the flux of SAX J1808.4-3658 between 2000 
March 5 and 8. With the XMM-Newton detec- 
tion of the source on 2000 March 6, we now have 
a clearer picture of how dim the source could be- 
come during certain phases of its 2000 outburst. 
As stated by Wijnands et al. (2001), the large lu- 
minosity variations observed during the 2000 out- 
burst are difficult to understand as due to sim- 
ilarly dramatic variations in the mass accretion 
rate. It is more likely that only modest variations 
in the accretion rate can trigger transitions be- 
tween two significantly different luminosity states. 
For example, centrifugal inhibition of accretion by 
the neutron star's magnetic field is expected be- 
low a certain critical accretion rate (the 'propeller 
regime'; Illarionov & Sunyaev 1975) and small 
but erratic variations in the accretion rate around 
this critical rate could in principle give rise to the 
enormous luminosity swings observed during the 
2000 outburst. When the source is in this pro- 
peller regime, accretion is inhibited, but it is evi- 
dent from the large brightness fluctuations during 
the 2000 outburst, that a considerable amount of 
matter was still available in the accretion disk. Ac- 
cording to Campana et al. (2002), a pure propeller 
contribution is ruled out in quiescence since this 
mechanism is expected to stop operating at lu- 
minosities below 10^^ erg s~^ because the source 
should turn on then as a radio pulsar. Stella et 
al. (2000) and Campana et al. (2002) suggested 
that a possible explanation for the quiescent flux 
is the emission from the shock front between the 
relativistic wind of the radio pulsar and the mat- 
ter out-flowing from the companion star. There- 
fore, the flux during the 2000 outburst observation 
could be higher than that observed in quiescence 
because of the large amount of matter still present 
close to the neutron star. 

Di Salvo & Burderi (2003) suggested an alter- 
native explanation for the quiescent emission of 
SAX J1808. 4-3658 in which this emission is pro- 
duced by direct dipole radiation from the radio 



pulsar. The quiescent X-ray luminosity and spec- 
tral shape of SAX J1808. 4-3658 are indeed con- 
sistent with those observed from several field mil- 
lisecond radio pulsars (see Becker & Pavlov 2002, 
for a recent review), but Grindlay et al. (2002) 
found that the millisecond pulsars in the globular 
cluster 47 Tuc have a predominantly soft spectral 
shape suggestive of a thermal origin in contrast 
to what is observed for SAX J1808.4-3658. Con- 
clusive proof for an active radio pulsar in SAX 
J1808. 4-3658 would come from the detection of ra- 
dio pulsations during quiescence, although such a 
detection might be inhibited by the ambient mat- 
ter still present in the system (see also the discus- 
sion in Di Salvo & Burderi 2003). 

Campana et al. (2002) noted that the quiescent 
properties of SAX J1808. 4-3658 are remarkably 
different from those observed for other quiescent 
neutron star systems: (a) its 0.5-10 keV luminos- 
ity is the lowest observed so far for any neutron 
star system, and (b) its spectrum is dominated by 
a power-law component instead of a thermal com- 
ponent. The quiescent emission of neutron star 
X-ray transients is most often explained by ther- 
mal emission from the neutron star surface releas- 
ing the heat deposited in the crust and core of 
the neutron star during outburst (see Brown, Bild- 
sten, & Rutledge 1998, and references therein). 
In this model, the exact luminosities of the sys- 
tems should depend on their time-averaged accre- 
tion rates (Campana et al. 1998; Brown, Bildsten, 
& Rutledge 1998). Due to the low peak lumi- 
nosity of SAX J1808.4-3658, Brown, Bildsten, & 
Rutledge (1998) predicted that this source should 
be rather faint in quiescence. At first sight the low 
detected quiescent luminosity is consistent with 
this prediction, but no strong evidence could be 
found for a thermal component in the quiescent 
spectrum obtained. This indicates that the ther- 
mal luminosity of this source is very low, implying 
a rapidly cooling neutron star which would require 
enhanced core cooling processes to explain (Cam- 
pana et al. 2002). 

The combination of the low time-averaged ac- 
cretion rate and the possibility of rapid core cool- 
ing might also be able to explain why the quiescent 
spectra of SAX J1808.4-3658 is dominated by the 
power-law component. If in the 'ordinary' quies- 
cent systems the thermal component were to drop 
to low luminosity levels similar to those of SAX 
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J1808. 4-3658, then their quiescent spectra would 
be dominated by the power-law component (as- 
suming that this component remains unchanged), 
just like in SAX J1808.4 3658. However, it re- 
mains to be determined if the power-law compo- 
nents in the different types of systems are due 
to the same mechanism, or if the quiescent spec- 
trum of SAX J1808.4-3658 is due to some other 
process, likely related to the different configura- 
tion and/or magnetic field strength of the neutron 
star in SAX J1808.4-3658. To get more insight 
into the nature of the unique quiescent proper- 
tics of SAX J1808. 4-3658, more neutron star tran- 
sients must be detected in quiescence. Particu- 
larly interesting systems are the two other, very re- 
cently discovered, accretion driven millisecond X- 
ray pulsars XTE J1751-305 (Markwardt & Swank 
2002) and XTE J0929-314 (Remillard, Swank, & 
Strohmayer 2002). It would be interesting to 
compare their quiescent properties with those of 
SAX J1808. 4-3658 and correlate possible similari- 
ties and/or differences in quiescence with the out- 
burst properties of those sources. 

This work was supported by NASA through 
Chandra Postdoctoral Fellowship grant number 
PF9-10010 awarded by CXC, which is operated 
by SAO for NASA under contract NAS8-39073. 
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Table 1 

Log of the XMM-Newton observation 



Observation ID 


Time of observation 


Instrumental modes" 


Live time'' 




(UTC, 6 March 2000) 




(ksec) 


0119940201 


18:29 - 23:53 


MOSl FF/Medium 


16.1 




17:09 - 20:24 


M0S2 SW/Medium 


11.4 


011!)!)4()r)()l 


21:o;-! - 2:-!:ol 


AI0S2 S\A'/McHlium 


:-!.!) 



*FF is Full Frame mode and SW is Small Window mode. Medium indicates 
that the medium filter was used during the observations. 

^Live time of the CCD on which SAX J1808.4-3658 is located, after eUmina- 
tion of the background flares. 
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Table 2 
Spectral parameters" 



Model 


(lO^i cm- 


Index/fcT 


Flux 

(10~i^ erg s~i cm~^) 


-^red 




Single 


component model 






Power law: 2000 
2001 

Black body: 2000 
2001 

Atmosphere'': 2000 
2001 


1 9+0.5 
'-■'^-0.7 

0.3+"-7 

< 0.1 
00+0.3 

< 1 
0.4+11 


2.2±0.3 

l-°-0.3 

0.37;H4 keV 
0.52t°;Ji keV 
0.20+|^ °^ keV 
0.27l°:°2 keV 


2.3 
0.62 
1.1 
0.29 
3.4 
0.24 


1.3 
1.1 

3.5 
3.2 
1.3 
4.4 


Multi component model, o 


nly for the 2000 outburst observation 




Black body + power law^ 

Atmosphere + power law'' 


1.3±0.5 
1±1 


0.2±0.1 keV 

1.6±0.3 
0.10l°:J^ keV 


2.6 
2.6 


1.2 
1.2 



'^The errors are for 90% confidence levels. The fluxes are unabsorbed and in the 0.5-10 
keV energy range. For the single component models, the fit parameters are given for the 
2000 outburst and also for the 2001 quiescent observations as reported by Campana et al. 
(2002). For the two-component models, only the 2000 outburst results are given. 

I'The atmosphere model by Zavlin, Pavlov, & Shibanov (1996) was used for the 2000 
outburst observations and with a distance of 2.5 kpc and a neutron star mass of 1.4 Mq. 
Campana et al. (2002) used the hydrogen atmosphere model by Gansicke, Braje, & Romani 
(2002) to obtain the fit result for the 2001 quiescent observation. 

'^The black body flux was 0.6 and the power law flux 2.0 x 10~i^ erg s~i cm~^ (0.5-10 
keV). 

'^The flux in the atmosphere component was 0.7 and the power law flux 1.9 x 10~i'^ erg 
s-i cm-2 (0.5-10 keV). 
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